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ABSTRACT 

We present deep optical integral-field spectroscopic observations of the nearby (z ~ 0.01) brightest 
cluster galaxy NGC 4696 in the core of the Centaurus Cluster, made with the Wide Field Spec- 
trograph (WiFeS) on the ANU 2.3m telescope at Siding Spring Observatory. We investigate the 
morphology, kinematics, and excitation of the emission-line filaments and discuss these in the context 
of a model of a minor merger. We suggest that the emission-line filaments in this object have their 
origin in the accretion of a gas-rich galaxy and that they are excited by v ~ 100 — 200 kms -1 shocks 
driven into the cool filament gas by the ram pressure of the transonic passage of the merging system 
through the hot halo gas of NGC 4696. 

Subject headings: Galaxies: elliptical and lenticular, cD — galaxies: individual (NGC 4696) — galaxies: 
interactions — galaxies:ISM — shock waves — techniques: imaging spectroscopy 



1. INTRODUCTION 

Many of the giant elliptical brightest cluster galaxies 
(BCGs) in the cores of massive clusters are surrounded by 
systems of filaments that produce a distinctive LINERQ- 
like emission spectrum. These structures exist in 'cool- 
core' clusters, which have centrally-peaked X-ray surface 
brightnesses. The densities and pressures are sufficiently 
high that the radiative-cooling timescales are less than 
a Hubble time in the core. As a result, there should 
be large-scale condensation and i nflow of cooling intra- 
cluster gas into the cluster center (fCowie fc Binnevlll977t 
iFabian fe Nulsenlll977t [Fabian et alJll984t l. However, re- 
cent observations show that much of the cooling of the 
cluster gas is being offs et by heating associated w ith the 
central galaxy (e.g., see iPeterson fe F abian 20061). 

In the current picture, the intra-cluster medium (ICM) 
cooling is regulated by intermittent energy feedback from 
the central galaxies. Potential heat sources include 
the mechanical energy of expanding radio lobes (e.g. 
iMcNamara fc N ulscn 20071 and r eferences therein) , con - 
duction fr om the cluster halo (e. g lVoigt fc Fabianll2004l ). 
hot stars llVoit fc Donahudll997f) . and high-energy parti- 
cles fe.g. lRephaeli fc SiHdll995f T Cool -core clusters allow 
detailed studies of galaxy and AGN feedback at relatively 
low redshift. They offer insight into analogous processes 
that occurred during galaxy growth in the high-redshift 
universe and into the physics needed to incorporate phys- 
ical descriptions of feedback in galaxy formation models. 

The extensive, filamentary emission-line nebulae in 
cool-core cluster BCGs are likely to be products of the 
ongoing and inter-related processes of heating and cool- 
ing. However, many questions about the origin of the 
material in the filaments and the mechanisms that ex- 
cite the emission spectrum remain largely unanswered. 

Various mechanisms have been proposed as sources 
of excitation in the filaments, including photoionization 
by radiation from the AGN, cluster X-rays, or hot stars 
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1 Low-ionization nuclear emission-line region (Hcckman 1980) 



(young or old); collisional heating by high-energy parti 
cles, shocks or cloud-cloud collisions, and conduction of 
heat fro m the X-ray coron a (see e.g. iFerland et a~ 
20091 : ICrawford fe Fabian) Il992t IDonahue et al 



200d lHatch et all 120071: Ijohnstone fc Fabianl 119881: 



Johnstone et al.ll2007t iSabra "et all 120001: iWilman et ail 
20021 and references therein). However, none so far 
satisfactorily describes the characteristic spectra, en- 
ergetics, and kinematics of the extended emission-line 
regions. 

Multi-wavelength observations of BCGs reveal that 
material in a range of phases co-exists in the filaments. 
Emission detected from various galaxies indicates the 
presence of cold CO at 10-100 K, mid-infrared-emitting 
H 2 at 300-400 K, molecular gas at 2000-3000 K producing 
near-infrared H 2 rovibrational lines, atomic gas at 10 K 
emitting atomic forbidden and hydrogen recombination 
lines, and a hot atomic com ponent traced by X-ray emis- 
sion at 10 6 - 10 7 K (see e.g. I Johnstone et al.|[2007l) . Ad- 
ditionally, these cooler gas phases are embedded in the 
hot, X-ray-emitting ICM at ~ 10 7 - 10 s K. Dust is 
also intermixed with the emission-line filaments in many 

lHatch et aiIl2007D . 

There are also systematic differences in the proper- 
ties of the nebulae amongst clusters. Surveys show 
that BCGs with more lumin ous line emission have: 
shorter cluster cooling times (10 'Pea et al.l 120081): in- 
creased blue /UV e xcess emission (iHeckman et al.l Il989t 
I Allen et al.l Il992t ICrawford et all 120051): greater IR 
luminosities (lEgami et al.l l2006at lO'Dea et all 120081: 
IQuillen et al.l 120081: IDonahue et al.l 120071 ) Har ger masses 
of m olecular gas associa ted with the filaments and galaxy 
(e.g. IP 'Pea et al.ll2008l): and more extensive emission- 
line regions jCrawford et al.l 119991 ) . Interestingly, there 
is also a clear trend of decreasing ratios of hydrogen re- 
combination to collisionally-excited forbidden line flux 
([N II], [S II]) with in creasing total emissio n-line lurni- 
edu.aft°llc«»^gnpatf.W (ICrawford et ai lE^p). Strong IR 
aim blue continuum excesses imply that significant star 



also intermix ed with the emission-lm e hlamcn ts m many 
systems (e.g. IDonahue fc Void 119931 ISparks et all 119891 : 
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formation occurs in some BCGs. However, in many cases 
the emission spectrum does not resemble that from a 
purely star-forming H II region. 

These results have led some authors to conclude 
that several excitation mechanisms may participate, 
and contribute t o vary ing degrees in different systems. 
iVoit fc Donahue! ([1997D suggested that sources of 'sup- 
plementary heating' produce the LINER-like proper- 
ties of the spectra, though not necessarily through the 
same mechanism in all systems. Similarly, Hatch et al.l 
(2006) propose that photoionization by hot young stars 
becomes increasingly dominant in high-luminosity sys- 
tems over an underlying excitation source that produces 
the low-ionization (LINER-like) spectrum and dominates 
the lower-luminosity nebulae. The nature of LINER- 
producing component(s) is uncertain, however. 

We present here a detailed study of the extended fila- 
ments surrounding the giant elliptical galaxy NGC 4696, 
the brightest cluster galaxy of the Centaurus cluster 
(Abell 3526). This is a relatively nearby system (z ~ 
0.01), in a cool-core clus ter with peaked X-r ay surface 
brightness in the center ([Mitchell et al.l[l975[ ). From a 
pure cooling-flow model the inferred mass deposition rate 
is -20 Moyr- 1 ([Fabian et al.l[l98l . 

Extend ed line emiss i on aro und NGC 4696 was iden- 
tified by iFabian et al.l (119821) . and the s t ructur es have 
been imaged in detail by iCr awford et al.l (f2005h . Their 
narrow-band Ha images reveal a complex system of fil- 
aments that extends over 50" (— lOkpc). The galaxy is 
a prototypical cool-core BCG in which the filaments can 
be spatially resolved and studied in some detail. 

The Ha luminosity of the galaxy emission places it at 
the low to intermediate end of the range of BCG line 
luminosity amo nst the ROSAT B rightest Cluster Sam- 
ple compiled by lAllen et al.l (|1992h . This sample spans 
four decades of Ha luminosity. The spectral properties 
are typical of lower-luminosity BCGs, exhibiting a very 
high ratio of [N II] to hydrogen recombination emission 
(John stone et all 119871: (Lewis et al.l 120031: iSparks et all 
11989( 1. Since in this case the low-ionization (LINER- 
like) properties dominate, we can investigate the nature 
of the LINER excitation with less contamination from 
other sources of excitation. 

The host galaxy is a giant elliptical central cluster 
galaxy, b ut an extended cD stellar halo has not been 
detected (|Mur et al.ll2006t IJerien fc Dresslerlll997D . The 
surface brightness profile is peaked within the central 
few arcseconds and there is evidence from the HST I- 
band image of dual n uclei within the c entral arcsecond 
of the core (~ 0.2kpc: lLaine et al.ll2003l) . The measured 
central velocity dispersion (257 kms -1 ) implies a cen- 
tral black hole mass on the order of (4 ± 1) x 10 s 
M (IRaffertv et al.ll200l using t he Magorrian Mbh — v 
relation ([Magorrian et all 119981: IGebhardt etTal! 120001) . 
NGC 4696 is associated with the low-power FR I radio 
source PKS 1246-41. The irregular lobes of radio emis- 
sion of this source appear to have been highly distort ed 
by the surrounding medium ([Taylor et al . 2006, 2 007D . 

Strong dust features are seen in the gala xy, including 
a conspicuous dust spiral (|Shobb rook 1966|) that can be 
traced from t he outer reg i ons o f the galaxy and in to the 
galaxy core (|Laine et al.l 12003ft . Close correspondence 
between the optical emission and dust features is ob- 



served in a number of BCGs, and this is well-illustrated 
in NGC 4696 where dust lanes are a ssociated with most 
of the extended opt ical filaments (Spa rks et al.l 119891 : 
ICrawford et al.l 120051 ) . The dust and line emission are 
concentrated to the south of the galaxy nucleus. The 
dust lane obscures the galaxy light, suggesting that it 
is positioned on the near side of the galaxy. That a 
large amount of dust exists in the filaments is a chal- 
lenge to models in which the cool gas condenses from the 
hot ICM, because dust is quickly, relative to the cool- 
ing timescale, dest royed by sputtering in the hot cluster 
halo environment (jVoit fc Donahuelll995l : [Donahue et al.l 

l2l)00h. 

ICrawford eTaTl ([200H) compare imaging of the ex- 
tended emission of NGC 4696 in the radio, optical, and 
X-ray wavebands. The complex system of optical line- 
emitting filaments is closely traced by soft (0.3 - 1 keV) 
X-ray emission. The striking resemblance between op- 
tical emission from warm gas at T — 10 4 K and X-rays 
from material at — 10 7 K is also seen in other BCGs, e.g. 
M87 and NGC 1275 ([Sparks et alJ l200i IFabian et al.l 
120031 ). and is anoth er challenge for models of the excita- 
tion. The images of lCrawford et all (2005) also illustrate 
that there is no clear correlation between the morphol- 
ogy of the radio emission and the optical filaments in this 
system. However, as observed in many cool-core clusters, 
the lobes of radio emission coincide with c avities in the 
high- energy X-ray emission from the ICM (|Fabian et al.1 
l2005h . 

In this paper we present the kinematics and excitation 
of the optical line emission detected in the filaments over 
the central 20" (~ 4kpc). Our data show material that 
is spiralling inward to the galaxy nucleus and uniform 
excitation properties across the emiss ion regions. The re - 
sults support the proposal made bv ISparks et al.l (1989) 
that the material in the emission-line filaments surround- 
ing this galaxy originated in a neighboring dwarf galaxy 
that has undergone a recent merger with the BCG. We 
describe a scenario in which the line emission is excited 
by low-velocity shocks that are driven by the motion of 
the infalling material. We present shock models that re- 
produce the emission-line flux ratios observed in the gas 
and surmise that a minor-merger is energetically capable 
of exciting the observed emission. 

Throughout this paper we adopt the cosmological 
parameters Hq = 71kms _1 Mpc -1 , flm — 0.27 and 
Q A = Q.73, based on the five-year WMAP results 
(IHinshaw et al.1 [2001). NGC 4696 h as a redshift of 
z = 0.0099 (|de Vaucouleurs et al.|[l991l ) and the spatial 
scale is 0.21kpcarcsec _1 . We assume a heliocentric dis- 
tance to the galaxy of 44Mpc. 

2. OBSERVATIONS 

The observations of NGC 4696 were obtained with 
the W ide-Field Spectrograph ('WiFeS: lDo"pita et al.ll2007l 
I2010D and the Australian National University (ANU) 
2.3m telescope at Siding Spring Observatory. WiFeS 
is a double-beam, image-slicing, integral-field spectro- 
graph that records optical spectra over a contiguous 
25 x 38" field-of-view. This spatial field is divided into 
twenty-five 1 "-wide 'slitlets', with 0.5" sampling along 
the 38" length. The position angle of the slitlets on the 
sky was 305° during the observations of NGC 4696. 

The WiFeS instrument is a dual-channel device with 
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separate gratings and cameras for the blue and red chan- 
nels. Each camera is equipped with 4kx4k pixel CCD de- 
tectors. In the low-resolution grating configuration that 
was used for these observations (employing the R3000 
and B3000 gratings and the RT560 dichroic), this pro- 
vides a continuous wavelength coverage from 3290 to 
9330 A in the reduced data. The spectral sampling is 
0.81 A/pixel in the blue data cube and 1.31 A/pixel in the 
red data cube. The average spectral resolution achieved 
was ~100 kms- 1 FWHM. 

The observations of NGC 4696 were made on 14-17 
April and 1-3 May 2009. The nights were clear and the 
seeing ranged between approximately l"and 3". Multi- 
ple 30 minute exposures were obtained for each data set, 
for total on-source integration times of 14 hours with 
the blue camera and 18.5 hours with the red (the blue 
detector was not available during part of the observing 
period). The pattern of science observations consisted of 
two pairs of galaxy frames separated by a 30 minute ex- 
posure of the sky. The scaled sky frame was subtracted 
from each of the four galaxy frames in this sequence. A 
dither pattern of 2" offsets was applied to the position of 
each galaxy frame. A nearby reference star was used to 
autoguide the galaxy exposures. 

Each set of four galaxy frames and one sky frame was 
preceded and followed by an observation of a spectropho- 
tometric standard star (LTT4364, HD128279, or EG131) 
and these were used to calibrate the absolute flux scale of 
the observations. Dwarf stars with spectra that resemble 
a featureless blackbody (including LTT4364 and EG131) 
were also measured during the night and used to remove 
atmospheric absorption features from the flux standard 
star and galaxy spectra. Arc calibration lamp frames and 
quartz lamp flat-field frames were observed at the begin- 
ning of each night, and twilight sky flat- field frames were 
obtained once during each of the two observing periods. 

2.1. Data reduction 

The data were r educed using the W iFeS IRAF data 
reduction package (Do pita et al.l I2010h ■ which is based 
on and utilizes the Gemini IRAF packagfl The pro- 
cedure began by treating each frame as a 2D image. 
The detector bias signals were removed in the standard 
way, using the overscan regions of the raw frames. Then 
the sky frame was scaled and subtracted from the bias- 
subtracted object frame. 

The spectrograph optically divides the field-of-view 
into slices along the long axis of the field. These 'slitlets' 
are stacked end-to-end and dispersed across the detector, 
to form a series of 2D slit spectra from adjacent slices on 
the sky. In processing, the 2D spectrum from each slice 
was extracted from the bias- and sky-subtracted frame 
and placed in an individual FITS image extension. Stan- 
dard methods for processing 2D spectra were then ap- 
plied individually within each slice as described below. 

The flat-field frame is formed from exposures of a 
quartz-iodine lamp. Within each slice, the spectra were 
normalized, using a spline function fit to the spectral pro- 
files. The shape of the lamp illumination in the spatial 
direction along the slices was corrected using twilight sky 
frames. An additional adjustment was made to maintain 
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the relative signal levels between the slices. The final 
flat-field frame then retained only pixel-to-pixel varia- 
tions in the spectral direction and the slit transmission 
profile in the spatial direction. These effects were sub- 
sequently removed from the science data by flat-fielding. 
Each slice of the galaxy frames was corrected with the 
corresponding 2D slice from the flat-field frame. 

The individual 2D spectra were then transformed to a 
common spatial and spectral coordinate grid. NeAr arc 
lamp frames were used to calibrate the spectral scales of 
the slices. The IRAF identify task was used to match 
the emission features in the arc lamp spectra, and fitco- 
ORDS was used to apply the spectral coordinate solution 
to the science data, slice by slice. The average RMS er- 
ror in the fit to the arc- line positions is 0.2 A for the red 
cube and 0.1 A for the blue cube. For the spatial align- 
ment, an occulting wire is available at an intermediate 
focus within the instrument. This provides a convenient 
spatial registration when the aperture is illuminated by 
the flat-field lamp. An image of the wire is seen near the 
center of each slitlet and this establishes the origin of the 
spatial coordinates in each slice. We assume a standard 
spatial scale of 0.5" per pixel along the slices. The data 
were resampled onto a common rectilinear grid by the 
transform task and the 2D slices were then stacked to- 
gether to form a 3D data cube. The ends of the spectra 
were trimmed from the reduced data cubes to remove the 
regions of low sensitivity in the grating+dichroic trans- 
mission profile. Approximately 1050 A and 550 A were 
removed from the short wavelength end of the red and 
blue cubes, respectively, and approximately 420 A and 
160 A from the long wavelength end of the cubes. This 
still provided ~ 500 A of overlap between the cubes. The 
trimmed regions do not extend to within 400 A of the 
emission lines measured. 

Telluric absorption features were corrected in the red 
galaxy spectra using an integrated spectrum from a re- 
duced smooth-spectrum standard star data cube. We 
normalized the spectrum to remove the continuum and 
obtain a spectrum containing only atmospheric absorp- 
tion features. This was compared to the absorption fea- 
tures seen in the galaxy spectra, scaled and shifted as 
required, and used to correct the galaxy data using the 
IRAF telluric task. After removing the telluric fea- 
tures from a standard star data cube, an integrated flux- 
standard spectrum was extracted and the standard IRAF 
flux-calibration tasks used to calibrate the spectra in the 
galaxy data cubes. The spectra across the whole field 
are corrected for the shape of the spectral transmission 
profile as measured at the position of the standard star 
(approximately the center of the field). The flux calibra- 
tion was performed without correcting for atmospheric 
extinction in the individual data frames. The variation 
in continuum flux among the frames gives an indication 
of the uncertainties in the flux calibration: there is a 
la variation of ~ 11%, attributable to the changing air- 
mass (the galaxy altitude was between 1 to 1.4 airmasses 
during the observations) and fluctuations in sky trans- 
parency. 

Each science frame was individually reduced by using 
the procedures described above. The dithered galaxy 
data cubes were then spatially registered using the con- 
soi-aHaHim centroid of the galaxy nucleus. This was mea- 



4 



Farage et al. 



sured for both the red and blue cubes, using an image 
formed by summing several spectral planes in a line-free 
region near the center of the spectrum. The cubes were 
then offset as required and median-combined with the 
IRAF task imcombine. Cosmic rays and detector de- 
fects were removed in this process. 

2.2. Continuum subtraction 

Measuring the flux of the extended line-emission re- 
gions surrounding NGC 4696 is complicated by the 
strong stellar continuum. This dominates the spectrum 
in most regions and must be removed to extract accurate 
emission-line fluxes. This can be done using the galaxy 
itself, rather than attempting to generate a theoretical 
stellar continuum template. 

Since the distribution of line emission around 
NGC 4696 is asymmetric (see Section 13.31) and the fil- 
aments do not extend beyond approximately 5 to 10" 
from the galaxy center to the north and east of the nu- 
cleus, we were able to extract a spectrum from a region 
of the inner galaxy that is free from line emission and use 
this to subtract the stellar component from the spectra 
at other locations. The representative galaxy spectrum is 
constructed from a 3"-diameter aperture placed ~ 8" to 
the north-east of the continuum center, where no line 
emission is detected in either the red or blue data cubes, 
as close to the galaxy core as possible. For each of the 
strong spectral emission features of interest, this spec- 
trum was scaled to match the flux level of the continuum 
in neighboring spectral regions (using the average level in 
two bands approximately 35 A wide bracketing the line 
or complex of lines and centered within 50 A of the fea- 
ture), before it was subtracted at each spatial pixel. An 
estimate of the error introduced in this process is made 
in the following section. 

2.3. Line profile fits 

After subtracting the galaxy stellar continuum, we ex- 
tracted the flux of the emission lines that are detected 
at better than the 3<r level. This corresponds to a line 
surface brightness >4x 10 -17 ergs -1 cm -2 arcsec -2 in 
the red data cube and >6x 10 -17 ergs -1 cm -2 arcsec -2 
in the blue for a line width of 300 kms -1 . The following 
lines are detected and fit: [O II] AA3726,29 A (spectrally 
unresolved), [O III] AA4959,5007 A, [N I] AA5198,200A 
(spectrally unresolved), [O I] 6300 A, [Nil] AA6548,83A, 
[S II] AA6716,31 A as well as Balmer hydrogen recombi- 
nation lines Ha, H/3 and H7. 

We fit single Gaussian line profiles to the emission lines 
of interest in the continuum-subtracted spectra using the 
IRAF task ngaussfit and measured line fluxes, the cen- 
troid positions, and widths at each spatial pixel having 
sufficient signal to allow this to be done. Ha and the 
lines of the [N II] AA6548,83 A doublet were fit together 
as three Gaussian components, and the lines of each 
of the [S II] AA6716,3lAand [O I] AA6300,63 A doublets 
were also fit simultaneously. All the other detected lines 
were fit individually. 

Radial velocities are given relative to the velocity of 
the galaxy with respect to the local standard of rest. The 
velocities have been corrected for the motion of the earth 
with respect to the local standard of rest; the magnitude 



of the correction required was less than 9.5 kms -1 for 
all observations. The instrumental resolution has been 
removed from the fitted line widths by subtracting in 
quadrature the average value of 97 kms -1 FWHM that 
was measured from an arc lamp calibration image. The 
results derived from the data cube and the emission-line 
fits are discussed in the following sections. 

Uncertainties in the relative line flux measurements 
were estimated using the noise properties of the spec- 
tra after continuum subtraction. For a given emission 
line, the la error in the flux density per pixel was taken 
to be the standard deviation in a line- free region adja- 
cent to the feature. These were combined in quadrature 
to form the uncertainties in the integrated line flux. 

3. RESULTS 
3.1. Galaxy continuum 

Figure [1] shows the integrated spectrum from a 5"- 
diameter circular aperture centered on the core of 
NGC 4696, prior to subtracting the galaxy continuum. 
Forbidden lines of N II, O II and S II, and the Ha recom- 
bination line are prominent, but the spectrum is domi- 
nated by a stellar continuum component that is typical 
of an early-type galaxy. The strength of the continuum 
increases toward longer wavelengths and exhibits absorp- 
tion features characteristic of evolved stellar populations: 
e.g. Ca K and H lines (at 3934 A and 3969 A, respec- 
tively), G band (CH+; 4300 A), blended Mg I b (5175 A), 
Na I D blended doublet (5893 A), and the broad TiO fea- 
ture (7126 A). In addition, the step in the average contin- 
uum intensity across the 4000A-break (.D4000 = 1-9) is 
typical of values in elliptical galaxies, including those of 
the low line-lumin osity BCGs in the RO SAT 'Brightest 
Cluster Sample' of ICrawford elaTl (fl99l . 

The galaxy spectrum shows no evidence for a blue 
stellar component that would indicate that significant 
star formation is occurring. This is in contrast to the 
blue and UV continuum excesses and large star forma- 
tion rates estimated in high-luminos ity BCGs: for exam- 
ple t hose in clusters Zw 3146 (e.g. IHicks k. M ushotzkv 
[2005]) . Abell 1835 (IMcNamara et all I2006D and Abell 
1068 ()Wise et al.ll2004[ ). It is consistent, though, with the 
general properties of the BGC population, which show 
redu ced signs of star formation in less luminous BCGs 
Ce.g.lAllen et al.lll992t[Crawford et aLlll999l:lEgami et al.l 
20063), such as NGC 4696 itself. 

We compare the radial surface brightness profile of the 
galaxy continuum with a Hernquist density model of the 
form 



following iLim et all (|2008|) . where M t is the total galaxy 
mass and a is a characteristic scale radius. We assume a 
constant /-band mass-to-light ratio of M/L = 4.7, taken 
from the correlation between stellar velocity di s persio n 
and /-band M/L presented by iCappellari et all ([20061 ) . 

A simple minimum least-squares criteria is used to 
compare the model with the measured galaxy surface 
brightness at 8000 A, averaged azimuthally between posi- 
tion angle 35° and 125° (i.e. in the quadrant of the WiFeS 
field where the dust absorption features are least promi- 
nent) and at radii from 2" to 17.5". The best fit is ob- 
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1. Integrated spectrum of NGC 4696 from a 5" diameter circular aperture centered on the galaxy nucleus, plotted against the rest 
gth. 
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Figure 2. Radial surface brightness profile of the stellar contin- 
uum emission m NGC 4696 at A = 8000A. Each point corresponds 
to the mean surface brightness of the pixels within a circular annu- 
lus. The error bars correspond to the dispersion in surface bright- 
ness of the pixels in each region. 



tained with a total mass M t = (4.0 ± 0.5) x 10 11 M© and 
scale radius a — 3.5 kpc. The measured and model radial 
s urface bright n ess pr ofiles are shown in Figure [5J 

lLaine et al.l ([2003) fit the /-band surface brightness 
profile with a 'Nuker profile' (a broken power law that 
has a turnover radius of 1.4") in the central 10", to fit 
the 'core' seen in the profile data. We find that the Hern- 
quist profile provides a better fit to the profile shape for 
our data beyond r ~ 10". 

Our estimated total stellar mass is different to the 
galaxy mass o f (11 .6 ± 0.1) x 10 11 M Q presented by 
Ra ffertv et al.l ([200(1 from i?-band imaging. The proce- 
dure used by the authors to determine this mass is not 
described in sufficient detail to understand the nature of 
this discrepancy. 

The stellar continuum is relativ ely redder close r to th e 
galaxy core, as also observed by Ide Jong et al.l (|1990f) . 
This may be either the effect of dust extinction increasing 
towards the galactic center, or else the result of a greater 
mean metallicity in the stars. 



3.2. Emission-line spectrum and line maps 

Figure [3] presents emission-line spectra in the region 
of the filaments where the line emission is strongest, be- 
fore and after continuum subtraction. The spectra are 
integrated over a 3"-diameter circular aperture centered 
approximately 3" north-west of the nucleus (AX = 0.2, 
Ay = 3.0 " in the reference frame of the images in Fig- 
ure |4j . The reference galaxy spectrum is shown, scaled 
to match the data in each panel, together with the total 
spectrum and the continuum-subtracted result, illustrat- 
ing the strongest emission-line detections. These spectra 
have not been corrected for reddening. 

Table Q] lists integrated, continuum-subtracted and 
extinction-corrected line-flux measurements in the same 
3" aperture, obtained from the Gaussian line profile fits. 
This table shows the emission features detected above 
the 3<T sensitivity level in the aperture. We correct the 
line flux measurement s for e xtinction using the redden- 
ing law of Osterbrock (1989]): the relative extinction at 
the wavelength of Ha is A-^a/Ay — 0.744 and at the 
wavelength of H/3 is Anp/Ay — 1.14. In addition, we 
adopt a value of R = Ay /E(B — V) = 3.1. To determine 
the total extinction, we assume that the intrinsic emit- 
ted flux ratio F (Ra) / F (R/3) = 3.1. This is somewhat 
higher than the classical Case B recombination value, 
but is appropriate for active galaxies, including LINERs 
(jVeilleux fc Osterbrock! [1987T ). in which there is a colli- 
sional excitat ion component in the Ha flux. From the 
dust maps of Schlcgel et al.l ([19981 ). the reddening con- 
tribution from our Galaxy in the direction of NGC 4696 
is A v ~ 0.34. 

The resulting extinction measured in the 3"-diameter 
aperture at the line peak is Ay = 0.74 ± 0.07 mag. The 
error is obtained from the uncertainty in the measured 
Ha/H/3 flux ratio integrated over the aperture. The mea- 
surements represented in Table Q] have been corrected for 
extinction of this magnitude. The measured Ha line flux 
in the aperture is (4.0 ± 0.4) x 10~ 15 ergs _1 cm -2 , 
and the extinction-corrected result is (6.7 ±0.9) x 10~ 15 
ergs -1 cm~ 2 . The error here is determined from the un- 
certainty in the flux calibration and the extinction. Line 
fluxes for the other detected emission lines are presented 
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Figure 3. These panels show regions of spectra near the strongest emission lines, extracted from a 3"-diameter aperture at the emission- 
line peak. The top (solid) line in each case is the raw integrated spectrum. Ovcrplotted in dotted lines are the line-free galaxy spectrum 
used to subtract the stellar component. The lower traces show the (integrated) emission-line spectrum produced by subtracting the galaxy 
continuum. 



as ratios to the Ha flux in Table [TJ Estimated uncertain- 
ties in the flux ratios, derived as described in Section [2751 
are also shown in the table. 

The top-right panel of Figure [3] indicates possible de- 
tections of the [Ne III] A3869A, H( A3889A and blended 
[Ne III] A3967A+ He A3970A lines, which are close to or 
below the 3<r detection threshold and could not be fit in 



Table 1 

Extinction-corrected flux ratios for emission lines integrated over 
a 3"-diameter aperture at the emission peak, relative to the 
measured Ha line flux: F(Ra) = (6.7 ±0.7) X 10 -15 
ergs -1 cm -2 . 



Emission line 
[O iTj 3726+3729 
H7 4340 
H/3 486f 

O III] 4959 

O III] 5007 

N I] 5198+5200 

O I] 6300 

N II] 6548 

N II] 6583 

S II] 6716 

S II] 6731 



F{X)/F(Ua) 
1.49 ± 0.06 
0.13 + 0.01 
0.32 + 0.01 
0.06 + 0.01 
0.31 + 0.01 
0.12 + 0.01 
0.20 + 0.01 
0.65 + 0.01 
2.05 ± 0.02 
0.61 + 0.01 
0.41 + 0.01 



individual pixels. 

Figure [4] presents maps of the emission- line properties 
derived from the Gaussian profile fits for three strong 
emission lines in the spectra: the line flux in the top 
panel (no reddening corrections have been applied), the 
velocity of the line centroid in the middle panel, and the 
derived velocity dispersion in the bottom panel. The 
cross marks the position of the galaxy nucleus, at the 
peak of the continuum emission surface brightness dis- 
tribution. 

3.3. Emission-line region morphology 

The distribution of the detected line emission from the 
filaments of NGC 4696 is illustrated in Figure [5j which 
shows a projection along the wavelength axis of the data 
cube over the [N II] AA6548,6583Aand Ha lines. The 
structures detected here can be matched to the feature s 
seen in narrow-band imaging bv lCrawford"et~aLl (120051) . 
The detected filaments extend to the east, south, and 
west of the galaxy center. Some filaments have a gen- 
erally radial orientation with respect to the galaxy, but 
there are curved arcs, bars, and filaments in the struc- 
tures that are extended at orientations that are much 
more nearly tangential. 



IFU Observations of NGC 4696 



7 



[N II] 65B3Aline flux 



Ha 6563A line flux 



[O III] 5Q07A line flux 




-5 5 10 
ax (arc sec.) 



-5 5 10 

ax (arc sec.) 



-5 5 10 

AX (arc sec.) 



Figure 4. Measured pixel flux, velocity centroid, and velocity dispersion of the [N II] A6583 A (left), Ha A6563 A (center), and [O III] 
A5007 A (right) lines. Note the strong increase in flux towards the center of the galaxy, the systematic rotation of the gaseous filaments, 
and the fairly uniform ~ 120 kms -1 velocity dispersion throughout the emission-line system. The apparent increase in velocity dispersion 
in regions of strong velocity gradient is probably due to beam-smearing in these areas. 



The most prominent feature is the bright filament 
that spirals from the nucleus to the north-west, and in 
the outer regions curves around the galaxy core to the 
east. This is the structure that was identified i n earl y 
imaging of the gas an d dust b y Fabia n et al.l (0 -982), 
iJorgensen et al.l (119831). jNorgaard-Nielsen fe; Jgrgensenl 
(|1984f ). and lSparks et all (|1989| ). The surface brightness 
of the filament decreases with distance from the galaxy 
nucleus along its length (other than at a brighter knot 
that appears at the intersection with another filament). 
The outer length of this main filament, to the south of the 
galaxy core, coincides closely with t he prominent dust 
lane that is seen c learly in the HST (lLaine et al.l I2003T ) 
and B — I color ()Crawford et al.l 120051) images of the 
galaxy. In the HST images, this lane of dust appears to 
be part of a continuous structure that spirals around the 
galaxy center at small radii to the north-east and with 



increasing radius to the west and south, passing through 
the brightest regions of the line emission. Lanes of ab- 
sorption also coincide with outer filaments to the south 
and west of the main filament. The association of the line 
emission and dust absorption features indicates that the 
observed system of gas and dust is inclined towards us, 
is on the near side of the galaxy, and has an intrinsically 
one-sided spiral distribution about the galaxy. 

The peak in surface brightness of the extended emis- 
sion occurs approximately 3" (0.6 kpc) to the north-west 
of the galaxy center. This peak occurs inside the curve 
of the dust lane. The spatial association of the main fila- 
ment with the dust lane is very clear in the outer extent 
of the filaments. Close to the nucleus, however, the main 
dust spiral departs from the ridge line of the line emis- 
sion along the filament. As seen in Figure [6l the filament 
of optical emission bends toward the nucleus from the 
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Figure 5. Maps of integrated N II 6548, 6583A+ Ho line emission from the continuum-subtracted NGC 4696 data cube. The image on 
the left shows the detail of the inner structures with a linear grayscale and contours of integrated flux. The right panel shows the same 
image with histogram-equalization of the levels to highlight the faint outer filaments. 



north-west, whereas the dust spiral continues around to 
the north of the core and finally turns towards the nu- 
cleus from the north-east. If the ridge of emission tracks 
a shock front, as proposed below, this might indicate that 
the shock front is progressively propagating through the 
dust lane as the dusty gas spirals in towards the nucleus. 

The morphology of the emission from weaker lines such 
as [O I] and [N I] is consistent with that of the strong 
emission lines (hence are not pictured in Figure [4]) . To 
the degree that our spatial resolution reveals, it appears 
that all the optical forbidden and hydrogen recombina- 
tion emission lines originate in the same gas. 

3.4. Line luminosity and mass of filaments 

The integral-field data permit us to determine the total 
emission-line luminosity of the filaments, which is useful 
to constrain both the mass and energy budgets. The flux 
measurements from the emission-line fits are summed 
over all the pixels where the value is greater than the 
3cr detection limit to obtain a total line flux. 

The total flux detected in Ha is (5.7 ± 0.6) x 10" 14 
ergs -1 cm -2 prior to correcting for extinction. The red- 
dening in the observations is estimated from measured 
Ha/H/3 line flux ratios as described in Section Since 
H/3 is relatively faint, we do not detect it over the same 
spatial extent as we measure the total Ha flux. This 
restricts our ability to map the spatial variation of the 
extinction. In spatial pixels where both lines are detected 
above the 5<r detection limit, the extinction correction is 
determined from the Ha /H/3 line ratio measured in that 
pixel. Everywhere else, we assume an average extinction 
level of Ay = 0.4 mag. This value is determined from 
the integrated line fluxes over the region where both Ha 
and H/3 are detected. The line luminosity is then derived 
by summing the individually extinction-corrected fluxes 
from each pixel. 

This results in an extinction-corrected, integrated Ha 
flux of (7.9±0.9) x 10 -14 ergs -1 cm -2 . At a distance of 



44Mpc, the Ha line luminosity is then (1.8 ±0.2) x 10 40 
ergs -1 . This is consistent w ith the lower-limit esti- 
mate of iCrawford et all ()2005l ) from narrow-band imag- 
ing (saturation in their images allowed only a lower 
limit to be obtained): L(Ra) > 1.5 x 10 40 ergs" 1 . 
The total [N II] A6583 A luminosity from our data is 
(3.7 ± 0.4) x 10 40 ergs -1 , si milar to the value o f ~ 
5 x 10 40 ergs -1 measured by Fab ian et al.l (|1982ft for 
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Figure 6. An HST image of the inner regions and dust lane in 
NGC 4696, obtained from the Hubble Legacy Archive (HST pro- 
posal 9427; dataset: HST_9427.06.ACS.WFC.F435W), with con- 
tours of optical emission ([N II] flux from our data) overlaid. The 
main optical filament is marked by the white points, which were 
selected by eye to follow the ridge of emission along the main fil- 
ament. This ridge line, initially along the outer edge of the dust 
lane, progressively moves inside the dust spiral towards the nucleus. 
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the [N II]A6583A flux. In earlier narrow-band work, 
Sparks et al. ( 1989) estimated an approximate combined 
Ha + [N II] A6548+6583A luminosity of 9 x 10 40 
ergs -1 , and given the inherent difficulties in calibrating 
such measurements and accurately extracting only the 
line emission, this is in reasonable agreement with our 
result. 

An estimate of the mass of emitting ionized hydro- 
gen can be made using th e Ha emission-line luminosity, 
L(Ha) (|Osterbrockll 19891) : 



Ah 



H+ 



L(Ha) m p 
n e a^ a h i/ H « 



(2) 



where we adopt the effective recombination coefficient 
for Ha line emission af^, = 1.17 x 10~ 13 cm 3 s _1 , from 
lOsterbrockl (| 1989( 1 for Case B emission at T = 10 4 K. 
Then 

M H+ ^4.3xW 100cm ^ M . (3) 



Thus, assuming an average electron density of ~ 
100 cm -3 , we find a total mass of ionized hydrogen in 
the filaments of M H + ~ (4.3 ± 0.5) x 10 5 M . 

The mass and density can be used to define a volume 
filling factor e for the ionized gas in the filaments. The 
volume filling factor for material in a structure is de- 
fined as the fraction of the encompassing volume that 
is filled by material at the estimated density. An as- 
sumption about the geometry of the emitting region is 
required. We can calculate a filling factor for the ionized 
material within a spherical volume that encompasses the 
filaments: 

„_M H+ /p 



i^ 3 



where M H + is the mass of ionized gas derived from the 
Ha flux, p is the mass density of the emitting gas, and 
R is the radius of a sphere encompassing the emitting 
region. Assuming an electron density n e ~ 100 cm~ 3 , 
R ~ 2.9 kpc and the ionized gas mass determined above, 
the volume filling factor is e ps 10~ 6 . This is compa- 
rable to values found for other filam ent systems (e.g. 
lOgrean et al.ll2010t lHatch et al.ll2007( l. and is consistent 
with the Ha-emitting gas being only a small fraction of 
the emitting material, located in a thin layer on the sur- 
face or within the volume of a larger mass of gas. 

This larger volume is unlikely to be spherical, so a more 
appropriate estimate might be obtained within a more 
constrained region of the filaments. We assume the emit- 
ting material from the 3"-diameter aperture used in Sec- 
tion [XU and Table [TJ at the line peak, is distributed over 
a spherical volume of radius R = 1.5" ~ 9.7 x 10 20 cm. 
From Equation [2] the estimated emitting mass of ionized 
gas in this aperture is M H + = (3.6 ± 0.5) x 10 4 and the 
volume filling factor is e « 10~ 4 , with n e ~ 100 cm~ 3 . 
This is still small, indicating a clumpy or filamentary 
distribution, but may more accurately reflect the volume 
fraction of the main filament that is occupied by ionized 
gas. 

The total gas mass is an important parameter, but is 
difficult to accurately ascertain. We discuss below var- 
ious estimates from the literature of the masses of the 
gas and dust in the filaments. These provide an indica- 



tionof the quantities of material present. iSparks et al.l 
(1989), for example, studied the dust absorption in op- 
tical imaging of the galaxy and calculate a mass of dust 
in the main dust lane > 2 x 10 7 M . This would then 
imply, using a Galactic gas to dust mass ratio, that the 
atomic component of the filaments has M~2x 10 9 M . 

Using the c olor g radient in the stellar continuum, 
Ide Jong et al.l (|1990l ) estimate the mass of dust present 
in the galaxy and dust lane from IRAS 100 /im observa- 
tions, by assuming that the gradient is entirely due to 
reddening rather than metallicity changes, and calculat- 
ing the mass of dust required to produce it. They obtain 
a total dust mass of ~ 0.2 — 5 x 10 6 M , with approxi- 
mately 10% of this mass located in the dust lane and the 
rest distributed over the central ^10 kpc of the galaxy. 
Assuming a Galactic gas-to-dust ratio, this implies a to- 
tal gas mass on the order of M ~ 3 x 10 8 M . 

Such estimates of dust mass cannot fully take into 
account the total gas mass, if the dust clouds become 
optically thick and therefore cool. The mass in BCG 
line-emitting filaments is believed to be dominat ed by 
molecular material in such cloud cores. lEdgel pOOlh have 
measured CO emission in a sample of luminous cool-core 
galaxies that is consistent with the presence of 10 9 — 
10 11 5 M of molecular gas at T ~ 40 K in these galaxies. 
Using Figure 10 of their paper, along with the inferred 
mass deposit i on rat e for NGC 4696 of 20 M yr _1 from 
iFabian et al.1 ()1982h then implies a molecular gas mass 
of M ~ 2 x 10 9 M in NGC 4696, in close agreement 
with the value inferred by ISparks et al.l (119891) from the 
dust obscuration. In NGC 4696. 1 Johnstone et all (|2007f l 
derive a mass of warm H2 of 1.3 x 10 5 M , from emis- 
sion within their Spitzer spectrograph aperture and esti- 
mate that approximately 10 times more mass may exist 
at cooler temperatures. 

As previously, we adopt an average reddening of Ay ~ 
0.4 along the extent of the main filament. Assuming 
a Galactic gas-to- dust ratio of N H /E (B — V) ~ 5.8 x 
10 21 cm" 2 mag" 1 (|Bohlin et al.lll978h . a length of 32", 
and width of 2" over which the emission in the filament 
is measured, we can make a rough estimate of the mass 
of material obscuring the line-emission region of M ~ 
10 7 M . This estimate is a lower limit for total mass 
because the ionized gas cannot be located entirely behind 
a dusty foreground screen, but is more likely mixed in 
with it. Over the total area of the line-emitting filaments, 
using the same average reddening, the resulting limit on 
the total mass is M > 8 x 10 7 M . 

Although these mass estimates are approximate and 
vary over a large range, they clearly indicate that there 
is a large mass of material in the extended filaments, of 
order a few times 10 9 M . It is most likely that this mass 
is largely in a cool, dusty, and dense molecular phase. 

3.5. Kinematics of the emitting gas 

There have been a number of long-slit spectroscopic 
studie s of the kinematics o f the gas in NGC 4 6 96. S tud- 
ies by ISparks et all (| 19971 ) and Ide Jong etaTl (| 19901 ) re- 
vealed smoothly varying line-of-sight velocities with val- 
ues chan ging by ^400 kms" 1 in the main filament. In 
addition, Spar ks et all (|1989l ) detect sodium absorption 
features in continuum-subtracted long-slit spectra that 
have similar kinematics to the emission-line gas along 
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both position angles observed. The measured line-of- 
sight velocities are between and 450 kms -1 . There is a 
cold, neutral component that is associated with the ion- 
ized filaments. They comment that this would not be 
expected in a cooling-flow model, but is consistent with 
a merger. 

Our integral-field spectroscopy greatly improves our 
ability to map the velocity dispersions and kinematics 
of the line-emitting gas. The middle panels of Figure [4] 
show the line-of-sight velocity of the emitting material 
measured from the Doppler shift of the Gaussian line pro- 
file centroids, relative to the galaxy. The bottom panel 
shows the line-of-sight velocity dispersion from the fit- 
ted profile widths, after correcting for the instrumental 
width (cr = 41 kms -1 ). 

The velocity dispersion detected across the main fil- 
ament is relatively uniform, with most of the material 
having velocity dispersions of ~ 120 ± 10 kms -1 (see 
Figure S]). Ridges of higher dispersion (160-200 kms -1 ) 
appear in the maps, but the dispersions in these regions 
are consistent with 'beam-smearing' caused by our fi- 
nite spatial resolution, which results in an increase in 
apparent velocity dispersion in regions of a steep spatial 
gradient in line-of-sight velocity. The measured velocity 
dispersions from the Ha and [N II] lines are effectively 
identical. 

The velocity dispersions in the main filament are sig- 
nificantly lower than the rotational/bulk velocities of the 
gas: without correcting for inclination, the ratio of rota- 
tion to turbulent velocity ranges from 1 to 2 in this struc- 
ture. The central stellar velocity dispersion in the galaxy 
is also a factor of two larger than the line-of-sight velocity 
dispersion of the line emission. The velocity dispersion is 
much larger than the thermal broadening (~7 kms -1 at 
10 4 K for hydrogen) and is therefore either generated by 
shocks or by the superposition of large-scale random tur- 
bulent motions. 

The main filament has broader lines and higher line- 
of-sight velocities than the outer filament arm to the 
south, in which the average velocity dispersion ranges 
from ^60 kms -1 to ^90 kms -1 and the line-of-sight ve- 
locity from ^100 kms -1 to ~150 kms -1 . This outer fil- 
ament also has lower surface brightness. 

The velocity distribution of the line-emitting gas shows 
an ordered structure with a component of rotation in the 
nebula. The velocity maps derived from the line centroid 
positions for other ions are essentially identical to the 
[N II] and Ha velocity maps shown in Figure^ indicating 
that the hydrogen recombination and forbidden lines in 
the optical spectrum originate from the same gas. 

The emission from the ionized gas is largely redshifted 
to the south of the galaxy center and blue-shifted to the 
north. There is a smooth variation in the line-of-sight 
velocity along the principal filament that traces the dust 
spiral. The velocity increases inward along the filament 
as it bends around the galaxy to the south of the core, 
and then decreases to become blue-shifted as the filament 
curves to the north past the galaxy nucleus. As a result, 
the kinematics of the gas in this filament resemble orbital 
motions: if the material is inf ailing, then the filament is 
inclined so that the outermost end is nearest to us and 
the material recedes along the filament to the south of 
the core, before turning around on the far side of the 



galaxy core onto an approaching trajectory to the north 
of the nucleus. 

At the outer end of the main filament at a radius of 
13" (2.7kpc), the average line-of-sight velocity of the 
emitting gas is 180 kms -1 . This velocity increases to 
280 kms -1 at a radius of 6.5" (1.4kpc) in the filament 
to the south-west. The maximum blueshifted velocity is 
150 kms -1 at a radius of 7" (1.5 kpc) to the north of the 
core. There is a total velocity range of 430 km s -1 in the 
me asured line c e ntroid s, similar to the value measured 
bv lSparks et~aTl (|1997| ). 

Assuming a Keplerian orbit with radius r = 1.5 kpc 
and circular velocity v — 300 kms -1 , the enclosed mass 
would be 3.1 x 10 10 Mq. The mass profile derived from 
the continuum surface brightness profile of the galaxy 
in Section 13.11 indicates that there is a stellar mass of 
3.6 x 10 10 M© inside this radius, in excellent agreement 
with our kinematic estimate, and showing that the con- 
tribution of dark matter to the total mass in this central 
region must be fairly small. The orbital timescale implied 
by material at a galactic radius of 1.5 kpc travelling at 
velocity ~ 300 kms -1 is ~ 3 x 10 7 yr. Taking the total 
gas mass (from Section l3~4|) as 2 x 10 9 M Q , the kinetic 
energy of the gas would be 1 /2Mv 2 ~ 1 x 10 57 erg. If this 
energy is radiated away by shocks on the timescale of a 
single orbit, the luminosity produced would be 2 x 10 42 
ergs -1 . The Ha line luminosity from our line measure- 
ments is (1.8 ± 0.2) x 10 40 ergs -1 , and from the shock 
models presented below, this is about 0.7% of the total 
luminosity. Thus, our estimated total shock luminosity 
is ~ 2.6 x 10 42 ergs -1 , in good agreement with our kine- 
matic estimate of this quantity. 

3.6. Spatial variation of excitation 

Previous long-slit spectroscopy of NGC 4696 has shown 
that the extended filamen ts produce emissio n with a high 
ratio of [N II] to Ha flux (|Lewis et al.ll2003D . We confirm 
this result. Figure [7] shows that the ratio of [N II] /Ha 
is roughly constant at 2.0 ± 0.1 with remarkably little 
variation in the value (~5%) across the inner filaments 
of the galaxy. This average value of the [N II] /Ha flux 
ratio and the Ha luminosity of the system reported in 
Section 13.41 are consistent with other BCG emiss ion sys- 
tems of similar luminosity ((Crawford et al.lll999[) . These 
results show that the characteristic line ratios such as 
the high [N II] flux relative to Ha and dominance of the 
other low-ionization species, are not simply nuclear phe- 
nomena confined to the galactic core, but persist over 
the extent of the detected emission in the filaments. In 
addition, the uniformity of the emission-line flux ratios 
across the inner filaments implies that whatever mecha- 
nism excites this spectrum is acting over the full extent 
(25" or 5 kpc) of the detected filaments. 

This finding contrasts with results from other BCG fil- 
ament systems. For example, lHatch et al.l (|2006l ) placed 
a spectrograph slit along several of the filaments in 
NGC 1275 and found the [N II] /Ha ratio decreases along 
the filaments with distance from the galaxy. The authors 
suggest that this may be the result of either decreasing 
metallicity or changing excitation, such as greater spe- 
cific star formation activity closer to the nucleus. Unlike 
NGC 4696, NGC 1275 is a very line-luminous system 
among BCGs and there is clear evidence of star-forming 
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cl ouds associated with t he extended filaments. 

Edward s et alj HOOS) compare optical IFU observa- 
tions of the extended emission in a sample of cooling-flow 
and non-cooling-flow BCGs. They find diverse emission 
properties and morphologies among their sample, includ- 
ing a range of behavior in the variation of the [N II] /Ha 
line flux ratios. In general, they attribute lower [N II] /Ha 
ratios to regions where young stars contribute most of the 
ionising flux. They conclude that the ionization mecha- 
nisms that produce the optical line emission are not uni- 
form amongst the sample and represent a complex inter- 
play between an AGN, cooling cluster gas, and star for- 
mation that may be triggered by AGN outbursts and/or 
tidal interactions with other galaxies. 

In IFU observations of a sample of other B CGs with 
a range of Ha luminosities, Hatch ct al. (2007) find dif- 
ferent degrees of [N II] /Ha variation in the observed sys- 
tems. For example, in RXJ 0821+0752 there is a dis- 
tinct correlation between the low [N II] /Ha and Ha flux, 
whereas in the BCG of Abell 2390 there is less varia- 
tion, but in general higher [N II] /Ha ratios are associ- 
ated with dust lanes and lower ratios with knots of blue 
emission. This may also reflect the trend amongst BCGs 
to display weaker forbidden emission in high-luminosity 
systems with abundant star forma tion (as evidenced by 
blue knots). iWihnan et al.1 (|2006l ) found little variation 
in the flux ratios over the measured extent of nebulae sur- 
rounding four high-luminosity systems. If, as suspected, 
such systems are dominated by more widespread regions 
of star formation, the emission-line ratios may be simply 
consistent with photoionization by hot stars. On this 
hypothesis, they suggested that in lower-luminosity sys- 
tems the flux ratios might show more spatial variation 
with higher [N II] /Ha ratios away from the less-abundant 
star-forming clumps. 

NGC 4696 may represent the low-luminosity extreme 
of this model, as there are no detectable star-forming re- 
gions associated with the filaments. In this case, the 
underlying low-ionization excitation mechanism domi- 
nates everywhere. M87, low line-luminosity BCG of 
the Virgo Cluster, may be a similar example, though 
in this case observations resolve the core of the galaxy 
and the emission-line ratios associated with the extended 
emission differ from those of the nucleus. The ex- 
tended fikinents_prc^uce_a_typically LIN ER-type spec- 
trum (jFord fc Butcherl Il979t ). found by iDopita et all 
(1997) to be consistent with shock models, whereas the 
higher-ionization emission from th e core is consiste nt 
with excitation by photoionization (Sa bra et alj |2003). 

The [O III]/H/3 ratio is a powerful diagnostic of the 
excitation of the gas. This ratio is also measurable in 
the central regions of the filaments in our data, though 
the signal-to-noise ratio is much reduced as a result of 
the overall weaker line emission. We find little evidence 
of a variation in the value over the inner filament. How- 
ever, there is evidence for a slight increase in the ratio 
at the galaxy core, where the value may be enhanced by 
radiation from the (weak) AGN in this system. This ef- 
fect can be seen in the line flux maps of Figure |4j as the 
peak of the [O III] emission is closer to the galaxy core 
than the peak of the hydrogen lines. However, over most 
of the extent of the filaments, the line ratios, including 
[O IH]/H/3, do not indicate a dominant contribution from 
AGN photoionization. The main filament uniformly pro- 



duces a low-ionization (LINER-like) emission spectrum 
(see Section [3~7l We do not detect He IA5876A, which is 
sometimes seen in more Ha-luminous BCG systems (e.g. 
iSabra et alJl2000f ). 

The intrinsic velocity dispersion in the ionized gas 
causes the [O II] A3729, 26 A doublet to be unresolved in 
our spectra. However, our wavelength range includes an- 
other density-sensitive pair of collisionally-excited lines: 
[S II] A6716,31 A. Figure maps the ratio of the flux in 
the 6731 A component to that in the 6716 A line. The av- 
erage value of the [S II] doublet flux ratio over the central 
filament is 1.5 ± 0.1, which is consistent with densities 
close to or within the low-density limit (< 300 cm -3 ). 
Again the [S II] flux ratio is uniform over the brightest 
regions in the filaments, with a la variation of 0.2 in the 
value. 

3.7. Diagnostic diagrams 

The most commonly-used optical diagnos t ic line -flux 
diagrams were developed by I Baldwin et al.l (|198lD and 
iVeilleux fc Osterbrockl(|1987h . These employ strong lines 
that are close in wavelength to decrease the effect of red- 
dening on the ratios. Such diagnostic diagrams are often 
used in classifying galaxies and identifying AGN versus 
starburst-dominated syste ms by distinguishing different 
ionization mechanisms (e.g. iKewlev et aLll2006f) . 

In Figure [3 we plot the flux ratios from Table [TJ at the 
emission-line peak on the common diagnostic diagrams. 
This figure shows the regions occupied by AGN, H II re- 
gi ons and LINERs, as defined by the classification scheme 
of IKewlev et al.l (I2006T) . developed from their earlier work 
(|Kewlev et al.ll2001[) . The typical line ratios from TableQ] 
for NGC 4696 are p lotted along with the i ntegrated line 
ratios measured bv iCrawford et all fl999) for a sample 
of X-ray-selected BCGs. The relative sizes of the points 
represent the relative Ha luminosities measured in that 
study. These demonstrate clearly the LINER-like prop- 
erties of the filament spectrum of NGC 4696, typical of 
a low Ha-luminosity system. Thus any conclusions that 
apply to the mode of excitation of NGC 4696 may be 
applicable to other low-luminosity BCGs. 

4. NGC 4696 AS A MINOR MERGER 

There is extensive work in the literature that addresses 
the question of how the mass of brightest cluster galaxies 
is assembled, including the fraction of the mass that may 
be obtained through merger and accretion events. The 
bulk of the mass of brightest cluster ellipticals is believed 
to have been assembled in this way. BCGs are located in 
the densest galactic environments and the galaxy lumi- 
nosities are weakly correlated with the clust er richness. 
X-ray luminosity, and velocity dispersion (e.g. lSchombertl 
119871 and references therein), signs that they are prod- 
ucts of their host cluster environment. There is also a 
high frequenc y of multiple nuclei amongst BCGs (e.g. 
IHoesselllniSOh . 

We have already accumulate d a variety of eviden ce to 
support the idea proposed by iSparks et al.l |l989) that 
the filaments of NGC 4696 represent the infalling rem- 
nants of a minor merger with a neighboring gas-rich clus- 
ter galaxy: 

• The kinematics of the emission-line gas show that 
the main spiral-shaped filament in NCC 4696 is a 
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Figure 7. Flux-ratio maps of the [N II] to the Ha recombination line emission (left), and for the density-sensitive [S II] doublet (right). 



coherent structure consisting of material with in- 
falling, orbital motions. 

The one-sided distribution of the emission-line 
streams is also readily explained by the accretion of 
another galaxy, infalling on a timescale comparable 
to the orbital timescale. 



• As ISparks et al.1 (|1989f) point out, the presence of 
large quantities of dust, and molecular, neutral, 
and ionized gas in the filaments is also consistent 
with a merger origin. 

• As discussed in Section 15741 there is evidence that 
the mass of material associated with the filaments 
is equivalent to that of a low-mass dwarf galaxy, 
of order a few times 10 9 M . However, this mass 
is much less than the inferred mass of NCC 4696, 
so the infall does not represent a major perturba- 
tion to the elliptical galaxy. It is therefore a minor 
merger. 

• Imaging of the core of NGC 4696 shows evidence 
for at least two nu clear components in the galaxy 
(|Laine et al.ll2003D . 

• The emission-line gas is characterized by a velocity 
dispersion of ~ 120 ± 10 kms -1 , characteristic of 
being shock-excited. 



• As discussed in Section 13751 the total emission- line 
luminosity from shocks is consistent with the esti- 
mated rate of dissipation of orbital energy in the 
filaments. 

In such a merger event, the dissipation of the orbital 
energy occurs through the drag forces generated by the 
motion of the merging galaxy gas through the hot galac- 
tic halo of the BCG. Shocks generated in the infalling 
gas would dissipate the orbital motion in this model, 
and therefore the signatures of shock excitation should 
be present. We investigate this possibility below. 



4.1. Radiative cloud shock models 

We have run a grid of radiative shock models using 
the photoionization and shock modelling code MAP- 
PINGS III, an updat e d vers ion of the code described in 
ISutherland fc Dopital (|1993D . The motivation for these 
simulations is a model in which the filaments are com- 
prised of clouds of cool gas that are remnants of a recently 
accreted, gas-rich, neighboring dwarf galaxy. These 
clouds are falling into the hot ISM halo of NGC 4696, 
which is compressed and heated by the ram pressure 
of the cloud's motion. The increased external pressure 
drives shocks into the cooler, dense cloud gas. The cool- 
ing time in the dense cloud material is short, so the cloud 
shocks are radiative and the observed optical line emis- 
sion is produced in the cooling region behind the shocks. 
Figure |H] shows a simple schematic of this model and is 
further discussed in Section l4~2l 

In the ID plane-parallel shock models, the pre-shock 
medium is taken to be the warm atomic component of the 
infalling gas, rather than the dense molecular medium 
that is visible in the dust absorption. The reason for this 
choice is as follows. The total surface brightness of a ra- 
diative shock of velocity v propagating into a medium of 
pre-shock density p scales as the mechanical energy flux 
through the shock S = 1/2/w 3 . For a system of shocks 
propagating into a fractal medium with a wide range of 
temperatures, the driving (ram) pressure P — pv 2 = con- 
stant. Thus v cx p -1 / 2 and therefore S oc p -1 / 2 . High 
density regions do not significantly contribute to the to- 
tal shock luminosity, not only because of this density- 
dependent factor, but also because the fractional filling 
factor of the denser cloud gas is small, reducing the total 
area that can radiate. Furthermore, in shocks propa- 
gating into molecular gas, the dominant luminosity is in 
collisionally-excited H2 infrared emission, not in atomic 
line emission. Indeed infrared emission lines from H2 gas 
are also seen, and the ir luminosity is correlated with the 
optic al line emission (jEdge et al.N2002t IJohnstone et al.l 
I2007T ). suggesting that these are produced by a common 
shock mechanism. 

The ionization state of the material prior to process- 
ing by the shocks (the pre-ionization) and the shock 
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Figure 8. The integrated line flux ratios from the 3"-diameter 
aperture located at the line peak in the main filament of NGC 4696 
are plotted here on optical diagnostic diagrams as filled rectangu- 
lar points, with side lengths that represent the uncertainties. For 
comparison, the flux ratios for ga laxies in ROSAT BCS clusters 
presented in Crawford ct al. (1999) arc also plotted. For the BCS 
points, the symbol size corresponds to the Ho luminosity of the 
system. 

velocity are treated as independent parameters in the 
models. The shock velocity is varied between 80 and 
220 kms -1 in intervals of 20 kms -1 , and the pre- 
ionization from neutral (i.e. the fraction of ionized hydro- 
gen is nnn/^H = 0) to fully ionized {nnu/nn — 1), with 
the fraction of hydrogen that is ionized in the pre-shock 
material increasing in steps of 0.1. A transverse mag- 
netic field that is consistent with equipartition of ther- 
mal and magnetic field energy is assumed (B — 5fiG for 
n = 10cm -3 ; B oc n -1 / 2 ). No external radiation field 



bow 
wave 




Figure 9. Geometry of the infalling cloud/minor- merger shock 
model. The infall velocity of the cloud and therefore the bow wave 
velocity is v\ and the cloud shock velocity is v s . Other labelled pa- 
rameters and components of the model are discussed in Section l4.2l 



Table 2 

Solar abundances Zq, the model abundance set log (n^/fi if), and 
included dust depletion factors log(D). 



Element 


Zq 


Models 


log(D) 


H 


1.00 


1.00 


0.00 


He 


-1.02 


-0.92 


0.00 


C 


-3.57 


-3.04 


-0.30 


N 


-4.17 


-3.34 


-0.20 





-3.31 


-3.01 


-0.23 


No 


-4.07 


-3.77 


0.00 


Ka 


-5.75 


-5.45 


-0.60 


Mg 


-4.43 


-4.13 


-1.15 


Al 


-5.56 


-5.36 


-1.44 


Si 


-4.49 


-4.19 


-0.88 


s 


-4.86 


-4.56 


-0.34 


CI 


-6.63 


-6.33 


-0.30 


Ar 


-5.60 


-5.30 


0.00 


Ca 


-5.69 


-5.39 


-2.52 


Fe 


-4.53 


-4.23 


-1.37 


Ni 


-5.79 


-5.49 


-1.40 



was added. 

The gas a bundances in Table [2] are taken from those 
presented bv lAsplund et al.l (|2009l ). with typical dust de- 
pletion factors, since dust physics is not treated explicitly 
in these models. Observations of the Centaurus clus- 
ter show that the ICM gas abundance in the vicinity of 
the galaxy and filaments is Z /Zq ~ 2 dGraham et al.l 
l2006USanders fe Fabianl 120061 ) and this overall metallic- 
ity is used in the models. However, we also find that 
an enhanced nitrogen abundance, by a factor of two, is 
required to better fit the observations. Nitrogen is made 
in the more massive AGB stars with hot-bottom burning 
(5 < M/Mq < 8), so possibly this nitrogen enhancement 
is associated with relatively less high-mass star formation 
and an older stellar population in the dense cluster envi- 
ronment. 

With this abundance set, the shock model is sought 
that provides a least-squares best fit to the values 
of the following line flux ratios: [O III] A5007/H,S; 
[N I] A5198+5200/H/3; [O I] A6300/Ha; [N II] A6583/Ha; 
and [S II] A6731/Ha. The individual line ratios from 
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Figure 10. Shock model results and measured line ratios for 
NGC 4696, plotted on the flux-ratio diagnostic diagram of [N II] 
6583A/Ha tw. [O III] 5007A/H/3. The solid lines connect the val- 
ues from models of the same shock velocity, from 80 kms -1 to 
220 kms -1 as labelled. The pre-ionization varies for the models 
that lie along these constant-velocity lines: the values correspond- 
ing to models with hydrogen pre-ionization fractions of 0, 0.3, 0.6, 
and 1.0 are marked by open circles that decrease in size with in- 
creasing ionization fraction. 



the grid of models are shown on the line flux ratio di- 
agrams in Figures [TUJ to [TSl along with the line ratios 
measured in the filaments of NGC 4696. The range of 
parameter space plotted here is much narrower than the 
normal range displayed on these type of diagnostic di- 
agrams. All except the [N I] A5198+5200/H/3 ratio are 
consistent with shock models with velocities in the range 
v s = 180 — 200km s _1 . We have determined the best-fit 
model spectrum by minimizing the RMS fit to all these 
ratios. The resulting best-fit model has a shock velocity 
v s = 180 kms -1 . The best-fit model also has complete 
pre-ionization of H, though it can be seen from the di- 
agnostic diagrams in Figures [TU] to [TS] that the results 
are not very sensitive to the pre-ionization of the high- 
velocity models. However, at shock velocities of 100- 
200 kms -1 the pre-shock gas will be photoionized by 
the upstream radiation flux from the shocked gas, so the 
most highly-ionized model is appropriate for the these 
velocities. The best-fit shock spectrum (at optical wave- 
lengths) is given in Table [3J 

Insufficient observed strength of the [O III] A4363A 
line has been used to argue against shock models as a pos- 
sible exciting mechanism in the filament s of BCGs (e.g. 
ISabra et al.|[2000t IVoit fc DonaTxul[l997l ). We do not de- 
tect the [O III] line in the spectra of NGC 4696, however 
this is consistent with the predictions of the shock mod- 
els described here. We can place an upper limit on the 
line strength of < 6 x 10~ 17 ergs -1 cm~ 2 arcsec -2 and 
therefore the ratio of fluxes [O III]/H/3 < 0.35, at the 
emission-line peak. From the model spectrum in Tabled 
the expected [O III] A4364Aflux is approximately 10% 
of the H/3 line strength, and so below the detectable level 
and consistent with the data. 

The only other emission line that is produced in the 
shock model spectrum at a level that may be detectable 
in our data is [Ne III] A3869 A; [Ne III] /H/3 = 0.34 from 
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Figure 11. Shock model results and measured line ratios for 
NGC 4696, plotted on the flux-ratio diagnostic diagram of [S II] 
6716+6731A/H a vs. [O III] 5007A/H^. Labels as in Figure [101 
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Figure 12. Shock model results and measured line ratios for 
NGC 4696, plotted on the flux-ratio diagnostic diagram of [O I] 
6300A/Ha vs. [O III] 5007A/H/3. Labels as in Figure llOl 



Table [3l However, as indicated above, at the blue end of 
the spectrum the ratio of detectable line flux to the flux 
in H/3 at the line-emission peak is > 0.35, so the predicted 
[Ne III] strength lies at the edge of our sensitivity limit. 
There is evidence for a marginal detection of this line 
at the emission peak (see Figure [3]), with an integrated 



flux of w (8 ±5) x 10~ 16 ergs 



(after extinction 



correction), but the fluxes in individual pixels are below 
the limit of detection. 

Though consistent in other measures, the shock models 
over-predict the strength of the [N I] A5200A doublet by 
a factor of ~1.4, the effect of which is seen in Figure [T3l 
Measurement errors may contribute to this discrepancy: 
the [N I] line flux is close to the detection limit in the 
blue spectra, and the line lies on the wing of the deep 
Mg I absorption feature, which could introduce system- 
atic errors in the galaxy continuum subtraction and line 
fitting that are not represented in the uncertainties given 
in Table [TJ Additionally, in the models the flux of the 
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Figure 13. Shock model results and measured line ratios for 
NGC 4696, plotted on the flux-ratio diagnostic diagram of [N I] 
5198+5200A/H/3 vs. [O III] 5007A/H/3. Labels as in Figure [TOl 




0.0 



0.2 0.4 0.6 

[OII]A3726,3729/ H/3 



0.8 



Figure 15. Shock model results and measured line ratios 
for NGC 4696, plotted on the flux-ratio diagnostic diagram of 
[O II] 3727A/HS vs. [N II] 6583A/Ha. Labels as in Figure [lOl 
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Figure 14. Shock model results and measured line ratios 
for NGC 4696, plotted on the flux-ratio diagnostic diagram of 
[O I] 6300A/Ha vs. [N II] 6583A/Ha. Labels as in Figure ITul 



[N I] line is very sensitive to the boundary conditions that 
terminate the simulations and so is strongly affected by 
the selected model configuration. 

In some sense, depending on projection effects, the ve- 
locity dispersion from the emission-line widths reflects 
the shock velocity. The line-flux ratios also depend on 
the shock velocity and should therefore correlate with 
the velocity dispersion. So if the line ratios are uniform, 
then the velocity dispersion is also expected to be uni- 
form. This is indeed observed to be the case in the main 
filament. 

We conclude then, that radiative shock models having 
super-solar abundances, and shock velocities of around 
v s = 180 kms -1 provide a good fit to the observed 
emission-line spectrum of NGC 4696 across the main fil- 
ament. We will now examine how we can use this cloud- 
shock model to further constrain the parameters of the 
infalling emission-line gas. 



Table 3 

The best-fit model shock spectrum 



A(A) 


F/F H0 


Ion 


A(A) 




Ion 


3425.81 


0.2341 


Ne V 


5015.70 


0.0286 


He 1 


3726.03 


2.0440 


O II 


5197.82 


0.2308 


N I 


3728.73 


2.4180 


O II 


5200.17 


0.2854 


N I 


3750.15 


0.0300 


H I 


5754.50 


0.0596 


N II 


3770.63 


0.0390 


H I 


5875.60 


0.1106 


He I 


3797.90 


0.0521 


H I 


6300.20 


0.3445 


O I 


3835.38 


0.0718 


H I 


6363.67 


0.1135 


O I 


3868.62 


0.3361 


Ne III 


6547.96 


1.8400 


N II 


3888.60 


0.0791 


He I 


6562.80 


2.9870 


H I 


3889.05 


0.1032 


H I 


6583.34 


5.4180 


N II 


3967.34 


0.1012 


Ne III 


6678.20 


0.0314 


He I 


3970.07 


0.1560 


H I 


6716.31 


1.3960 


S II 


4026.20 


0.1757 


He I 


6730.68 


1.1790 


S II 


4068.50 


0.0694 


S II 


7065.20 


0.0111 


He I 


4076.27 


0.0229 


S II 


7135.67 


0.0571 


Ar III 


4101.73 


0.2539 


H I 


7319.65 


0.0201 


O II 


4340.46 


0.4614 


H I 


7320.24 


0.0620 


O II 


4363.15 


0.0983 


O III 


7330.10 


0.0333 


O II 


4471.50 


0.0378 


He I 


7330.70 


0.0325 


O II 


4566.78 


0.1473 


Mg I 


7750.95 


0.0137 


Ar III 


4685.74 


0.0628 


He II 


8617.00 


0.0142 


Fe II 


4714.10 


0.0133 


Ne IV 


8727.66 


0.0149 


C I 


4724.11 


0.0144 


Ne IV 


8862.79 


0.0142 


H I 


4725.60 


0.0133 


Ne IV 


9014.91 


0.0189 


H I 


4861.32 


1.0000 


H I 


9068.44 


0.0956 


S III 


4958.83 


0.3342 


O III 


9229.02 


0.0262 


H I 


5006.77 


0.9627 


O III 


9530.44 


0.2370 


S III 



4.2. Global model parameters 

In our model, the infalling motion of the gas clouds 
generates the ram pressure that drives radiative shocks 
into the cloud material. The cloud velocity and the prop- 
erties of the hot, ionized gas in the galaxy halo (ISM) 
determine the magnitude of the driving pressure that is 
obtained. 

Figure [9] outlines the geometry of the model and the 
components referred to in the following discussion. As a 
cool gas cloud falls through the ambient medium, a bow 
wave, or shock if supersonic, forms in the ISM ahead of 
the cloud. The properties of the ambient ISM gas are 
denoted with the subscript 1, and v\ is the velocity of 
the ISM relative to the cloud (i.e. having magnitude 
equal to the infall velocity). Immediately behind the 
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bow wave the gas has pressure, density, and temperature 
Pi-, Pi, T2, and the velocity relative to the bow wave is 
v 2 . If v\ is greater than the sound speed in the ISM 
c S; i = (7Pi/pi) 1 ^ 2 , the bow wave is a shock, and the 
Mach number is Mi = «i/c Sj i > 1. In this case the pres- 
sure will increase across the shock front. If the motion 
is subsonic V\ < c S) i, so M\ < 1, and the pressure is 
constant across the bow wave p 2 = Pi- 

We determine the maximum pressure that acts on the 
cloud (due to its motion through the ISM) as the stag- 
nation pressure at the leading surface of the cloud where 
the flow comes to rest. The Bernoulli equation along a 
streamline traversing the post-bow-wave region, assum- 
ing an adiabatic, compressible flow behind the bow wave, 
gives: 



1 2 



hv = Streamline Constant = — n 



where v s — is the velocity at the stagnation point, and 
h is the specific enthalpy, 

, 7 P c 2 s 



7 — 1 p 7 — 1 
From this, the stagnation pressure can be expressed as 



Ps = Pi 



1 



7-1 



Mo 



7/(7-1) 



(5) 



in terms of the conditions in the flow behind the bow 
wave. We compare the infall velocity of the cloud to the 
ambient sound speed in the ISM to estimate the Mach 
number of the bow wave (Mi) and hence the properties 
o f the gas behind it. 

I Graham et al.l (|2006[ ) present density and temperature 
profiles of the IGM gas derived from X-ray observations 
obtained with Chandra, XMM- Newton, and ROSAT. We 
use these to estimate the pressure and sound speed in the 
halo of NGC 4696 and constrain the pressure driving the 
shocks in the filament clouds. At a radius of ~ 1 kpc (5"), 
which corresponds to the inner region of the filaments, 
the X-ray-emitting plasma has a temperature of T\ = 
(8 ± 3) x 10 6 K and hydrogen particle density Uh,x = 
0.13 ± 0.04 cm -3 . The sound speed in this gas is then 
c s .i = 440±50 kms" 1 and the pressure is pi/k — (2.5± 
0.5) x 10 6 K cm -3 . 

From our IFU observations of the inner filament, we 
assume an infall velocity of the clouds in the filament of 
vi ~ 300 kms -1 , approximately the maximum line-of- 
sight velocity of the emission in the filament. The Mach 
number of the flow is then Mi ~ 0.7. In this case the flow 
is in the transonic regime and a weak bow shock forms. 
For a transonic shock Mi « 1 w M 2 , and the pressure 
remains approximately constant across the shock: pi w 
p 2 ~ piVi/^f. From Equation [SJ the expression for p s 
becomes 



P1V1 

7 



7 + 1 



7/(7-1) 



: 1.23 pivl = 2.05 pi for 



7= 3' 



(6) 
(7) 



so the static pressure of the X-ray-emitting ISM is en- 
hanced by a factor of approximately two at the surface 
of the cloud. 



At the radius of the inner filaments, the mass density of 
the ISM is pi = lAm H n HA ~ 3.0 x 10 -25 cm -3 . With 
a bow-shock velocity of vi ~ 300 kms -1 , the resulting 
stagnation pressure is p s /k = 2.4 x 10 6 K cm -3 . This 
post-bow-shock pressure drives the cloud shock by pro- 
viding an equivalent pressure to the gas behind the cloud 
shock. We assume this pressure remains approximately 
constant across the post-cloud-shock region of the cloud 
(denoted by subscript 3 in Figure [9]). 

For a strong shock, the pressure in the post-cloud- 
shock region p% is related to the pre-shock cloud gas 
properties (p c , p c ) as follows: 



P3 : 



7 + 1 
2 

7 + 1 



PcV s 



P m H n H ^ c v s 



(8) 
(9) 



and the pre-shock hydrogen particle density in the cloud 
is then 



, 7 + 1 Ps 
2 p ran v 2 

a. 76 x 10 9 p s 



"//.,' : -— : T (10) 

(11) 

for 7 = 5/3 and the best- fit shock- model velocity v s — 
180 kms -1 . Hydrogen is fully ionized in the pre-shock 
cloud gas, so p — 1.4. From these results we can estimate 
the density of the ionized cloud material. For the inner 
filaments we then find a hydrogen density nn, c ~ 1 cm -3 . 

The pre-shock density in these models cannot be con- 
strained by the [S II] A6716/6731 Aratio, as the obser- 
vations are consistent with the low density limit (n e < 
300 cm -3 in the recombination zone of the shock, with 
T e ~ 6000 K). This prevents a direct measurement of the 
ram pressure associated with the shock. A somewhat less 
accurate estimate can be obtained from the Ha surface 
brightness, since this scales as the total emitted shock 
surface luminosity. The best-fit shock model has emit- 
ted surface luminosity in Ha: 



Sua 



"total 



-pv 



= 9.6 x 10" 



ergs 1 cm 2 . 



(12) 

where n e ^ c is the electron density in the pre-shock cloud 
material. 

The average Ha surface brightness in the main fil- 
ament measured in the 3"-diameter aperture around 
the region of the line peak is I^a ~ 9.4 x 
10 -16 ergs -1 cm -2 arcsec -2 . The measured surface 
brightness varies across the filaments, with a typical 
value in the outer region of the main filament of ~ 
1 x 10 -16 ergs -1 cm -2 arcsec -2 . This corresponds to a 
measured surface luminosity from the inner filament of 



'He 



I 



AttD 2 ( 206264. 



Ho 



V 



D 



5.35 x 10 u / Hct 



Q 3 ) 

where <p is the filling fraction of shock emission in- 
side the 3" aperture. This fraction is highly uncer- 
tain. However, if we assume 0=1, we obtain Sn a ~ 
5.2 x 10 -4 ergs -1 cm -2 in the central filaments and 
~ 5.5 x 10 -5 ergs -1 cm -2 in the outer region. Us- 
ing the relation between the density and surface lu- 
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minosity in Equation I12[ we find a pre-shock density, 
n e ,c ~ 5 cm -3 near the nucleus, and n e>c ~ 0.6 cm~ 3 in 
the outer regions. In the ionized pre-shock medium, 
the equivalent hydrogen densities are ~ 4 cm -3 and 
~ 0.5 cm~ 3 . There are uncertainties associated with 
these estimates, most notably in the filling factor of the 
shock surfaces in the surface luminosity estimates. How- 
ever, the result is consistent with the density estimates 
of ~ lcm~ 3 derived through the bow-shock/cloud-shock 
pressure arguments presented previously. The radia- 
tive shock model appropriate to the nuclear region with 
v s = 180 kins -1 , no ~ 4 cm -3 , and equipartition mag- 
netic field B — 5 fiG produces a mean density in the 
[S II]-emitting zone of the shock of 190 cm -3 , consistent 
with the [S II] A6716A/673lA flux ratio being observed 
at its low-density limit. These results provide support 
for the feasibility of the model as a mechanism for gener- 
ating shock-excited emission-line spectra in the extended 
filaments of NGC 4696. 

The increased pressure in the bow shock region will 
lead to an enhancement of the X-ray emissivity in 
the vicinity of filaments. This can provide a natu- 
ral explanation of the correlation bet ween the X-ray 
and t he optical emission observed by iCrawford et al.l 
(2005). These authors measure the 0.1 — lOkcV flux 
in the filament 22" to the south-east of the nucleus to 
be ~ 7 x 10 -15 ergs -1 cm -2 arcsec -2 . From our data, 
the Ha surface brightness in this outer filament is ~ 
5 x 10 -17 erg s - 1 cm -2 arcsec -2 , so the Ha flux in the 
region is approximately 0.7% of the soft X-ray flux. The 
total surface brightness in the cloud shocks, from the 
model, is about 150 times the Ha surface brightness. 
Combining this with the measured ratio of X-ray to Ha 
flux, the ratio of the X-ray to total cloud shock surface 
brightness is estimated to be Sx/S c ~ 0.9. Thus the 
energy dissipation in the X-ray plasma is comparable to 
that of the cloud shocks. This ratio might be lower if 
there are faster (non-radiative) shocks in the infalling 
gas, as these could provide soft X-ray luminosity associ- 
ated with the infalling gas, rather than the shocked ISM. 
Nonetheless, it is clear that the gas-dynamic drag of the 
hot ISM is an important factor in determining the accre- 
tion timescale of the infalling gas towards the nucleus. 

4.3. Extended LINER emission and shocks in other 
galactic environments 

Although LINER emission in galaxies is frequently a 
nuclear phenomenon, there are a range of galactic envi- 
ronments that host extended, extra-nuclear, LINER-like 
emission regions. Other than cool-core cluster BCGs, 
these include galactic winds, AGN jets, and interacting 
or merging galaxies. 

A link between LINER-like emission and shock exci- 
tation was proposed in early studies of LINER galaxies 
and has since been extensively di scussed and debated 
in the literature (e.g. iDopital 119941 ) . Indeed, from early 
invest i gations of LINER sp ectra . IKoski fe Osterbrockl 
(|1976f) . iFosburv et al.l (| 19781 ). and IHeckmanl (|1980D corT 
eluded that shock heating is responsible for producing 
the low-ion i zation spectrum in LINERs. More recently, 
iFiUppenkol (|2003T ) noted that despite a lack of strong ev- 
idence for shock heating in galactic nuclei, shocks are 
likely to dominate in cases of spatially-extended LINER- 



like ionization. A connection between LINER emission 
and shock excitation has been propo sed in models o f 
galaxy 'superwi nds'. e.g. in NGC 839 jRichet^ [2010l) 
and NGC 3079 (IFilippenko fc Sargentill992l). and i n in- 
teractions and mergers (e.g. IDelgado fc Perezl [19961) . 

Most interesting in the context of this paper, are recent 
findings that increasingly point to shock excitation of 
LINER-like emission i n the tidal r egions of galaxy merg- 
ers and interactions. lOgle et ail ()2007t ) present Spitzer 
observations of the luminous infrared emission associ- 
ated with the double galaxy system 3C 326. The north- 
ern galaxy of the pair produces strong H2 line emission 
and a LINER-type optical spectrum. The authors in- 
vestigate several heating mechanisms for the IR emis- 
sion and conclude that the molecular hydrogen spectrum 
may be powered by accretion shocks that arise in inflows 
induced by tidal int e ractions bet ween the two galaxies. 
IMonreal-Ibero et al.l (|2010l I2006D investigate the ioniza- 
tion of the extended emission regions in a sample of lumi- 
nous infrared galaxies (LIRGs) using integral-field spec- 
troscopy. They find that the emission in isolated systems 
is generally consistent with ionization by radiation from 
young stars. However, large fractions of the emission 
spectra from the extended regions in interacting pairs 
and more advanced mergers are consistent with ioniza- 
tion by shocks with velocities of 150 — 500 kms" 1 . In 
their model, the shocks are generated by tidally-induced, 
high-velocity flo ws that can be a resu lt of the merg- 
ing process (e.g. IMcDowell etaLl 12001 . They suggest 
that tidal forces may be central to the origin of shocks 
that ionize the gas in the extended regions and pro- 
duce the low-ionization spectrum. The most luminous 
BCGs do produce IR luminosities t hat place them in th e 
class of LIRGs (L m > 10 n L Q ; e.g. lEgami et~aLll2006bj) . 
While NGC 4696 is not such an IR-lumi nous system: 
L m = 2.3 x 1O 9 L (jKaneda et al J 120071) . results like 
those discussed above support the suggestion that infall 
and merger activity provides a mechanism for generating 
galaxy-wide LINER-like emission through intermediate- 
velocity shock excitation. 

Lastly, we note the observations presented by 
lEdwards et~aT1 (|2009f ) of the BCG system in the Ophi- 
uchus cluster, which make an interesting comparison 
with our results for NGC 4696. In Ophiuchus, the BCG 
does not produce detectable line emission, but is accom- 
panied by a nearby extended object that emits strong 
low-ionization optical lines. This object is ~ 1.7kpc from 
the center of the BCG on the sky and is separated by 
600 kms -1 in radial velocity. The red emission line ra- 
tios fall in the AGN or LINER regions of the diagnostic 
diagrams, but the [O HI]/H/3 ratio is not available to 
distinguish between these classes. The authors conclude 
that it is most likely that the emission is produced within 
a large gas cloud falling towards the BCG through the 
galaxy halo. It is possible that the emission in this struc- 
ture is excited by the ram pressure of its motion through 
the hot atmosphere of the BCG, as in the model we have 
discussed here. The Ophiuchus system may be in the 
early stages of a process of minor merger like the one 
that has occurred in NGC 4696. Further study is needed 
to better understand the nature of the interaction. 

5. CONCLUSION 
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Our IFU observation of the extended emission sur- 
rounding NGC 4696 provides a picture of material falling 
in along the bright curved filament that encircles the core 
of the galaxy, with an emission- line spectrum - and hence 
excitation mechanism - that exhibits uniform properties 
over the inner filaments. 

The kinematics indicate 'orbital' trajectories, spiralling 
in toward the galaxy nucleus. In addition to other obser- 
vational data that provide evidence for a recent merger 
event in the system, this prompts us to investigate a 
model in which the extended emission-line filaments form 
from the stripping of material during a minor merger 
event that also provides the source of energy that excites 
the line emission. In this model, the infalling clouds in 
the filaments are remnants of a low-mass, gas-rich (con- 
taining approximately 10 9 M a gas) cluster dwarf galaxy 
captured in the potential of the BCG. The motion of 
the filaments through the hot ICM produces bow waves 
that drive shocks into the clouds themselves. We present 
MAPPINGS models of the cooling behind these shocks 
and find that the optical emission-line spectrum is well- 
matched by ^180 kins -1 shocks propagating into ion- 
ized gas. 

The minor-merger scenario described here has the ben- 
efit of naturally accounting for the high densities and 
pressures in the filaments relative to the surrounding 
medium (unlike direct cooling models that assume that 
the filaments are in pressure equilibrium with the envi- 
ronment). It also explains the kinematic properties of 
the filaments, and is energetically self-consistent. Fur- 
thermore, the presence of large quantities of dust closely 
associated with the optical filaments is consistent with 
the model of an infalling satellite. 

Excitation by low- velocity shocks is a mechanism that 
may be capable of explaining the underlying LINER-likc 
properties of the extended filament emission in central 
cluster galaxies. Although such shocks may not be al- 
ways driven by infall or mergers, there is an interest- 
ing parallel between the properties of the filaments in 
NGC 4696 and those of emission regions in LIRGS, which 
have been attributed to tidally-induced accretion shocks. 
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